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Abstract 

We study three dimensional oceanic Lagrangian Coherent Structures (LCSs) in the Benguela region, as obtained from an output 
of the ROMS model. To do that we first compute Finite-Size Lyapunov exponent (FSLE) fields in the region volume, characterizing 
mesoscale stirring and mixing. Average FSLE values show a general decreasing trend with depth, but there is a local maximum at 
about 100 m depth. LCSs are extracted as ridges of the calculated FSLE fields. They present a "curtain-like" geometry in which 
^ the strongest attracting and repelling structures appear as quasivertical surfaces. LCSs around a particular cyclonic eddy, pinched 
off' from the upwelling front are also calculated. The LCSs are confirmed to provide pathways and barriers to transport in and out 
GN of the eddy. 

Keywords: Lagrangian Coherent Structures, Finite-Size Lyapunov exponents, ocean transport, Benguela upwelling region, 
oceanic eddy 

o 

^ ^ Mixing and transport processes are fundamental to determine 
^ the physical, chemical and biological properties of the oceans. 



1. Introduction. 



^^From plankton dynamics to the evolution of pollutant spills, 
^ ^ there is a wide range of practical issues that benefit from a 
correct understanding and modeling of these processes. Al- 
though mixing and transport in the oceans occur in a wide range 
of scales, mesoscale and sub-mesoscale variability are known 



o play a very important role ( Thomas et al. 2008 Klein and 





Q\ |Lapeyre|[20Q9 ). 
r>" Mesoscale eddies are especially important in this aspect be- 
cause of their long life in oceanic ffows, and their stirring and 
mixing properties. In the southern Benguela, for instance, cy- 
clonic eddies shed from the Agulhas current can transport and ^° 
exchange warm waters from the Indian Ocean to the South 
Atlantic ( [Byrne et alT| |1995t [Lehahn et aT) |2011i ). Moreover, 
mesoscale eddies have been shown to drive important biogeo- 
chemical processes in the ocean such as the vertical ffux of 
nutrients into the euphotic zone ( [McGillicuddy et aL) |1998[ 45 
^ [Oschlies and Gargon] |1998| ). Another eff'ect of eddy activity 46 
seems to be the intensification of mesoscale and sub-mesoscale 47 
variability due to the filamentation process where strong tracer 48 
gradients are created by the stretching of tracers in the shear- 49 
ipnd st rain-dominated regions in between eddy cores fElhm aidi] 50 
[et al.[ |1993| ). Studies of the vertical structure of such eddies 51 



et al.| ( [2009| )) have shown that they can extended to one thou- 
sand meters deep waters. 

In the last decades new developments in the description and 
modelling of oceanic mixing and transport from a Lagrangian 
viewpoint have emerged ( [Mariano et al.[|2002[|Lacasce 2008| ). 
These Lagrangian approaches have become more and more fre- 
quent due to the increased availability of detailed knowledge of 
the velocity field from Lagrangian drifters, satellite measure- 
ments and computer models. In particular, the very relevant 
concept of Lagrangian Coherent Structure (LCS) ( |Haller[[2000[ 



Haller and Yuan 2000) is becoming crucial for the analysis of 



• 16. 

X 



transport in ffows. LCSs are structures that separate regions of 
the ffow with diflferent dynamical behavior. They give a general 
geometric view of the dynamics, acting as a (time-dependent) 
roadmap for the ffow. They are templates serving as proxies to, 
for instance, barriers and avenues to transport or eddy bound- 
aries ([Boff'etta et al. 2001; Haller and Yuan 2000; Halierj [2002} 
[d'Ovidio et al. ^.2004,, 2009.; Mancho et al.|,2006) . 

The relevance of the three-dimensional structure of LCSs be- 



ins to be unveiled in atmospheric contexts ( du Toit and Mars 



26 jn the Benguela region (e. g. [Doglioli et al.| ( [2007] ) and [Rubio 



pen] [20Tq1 [Tang et al.|[20TT|[Tal lapragada et al.| |2011| ). In the 

case of oceanic ffows, however, the identification of the LCSs 
and the study of their role on biogeochemical tracers transport 
has been mostly restricted to the marine surface ( d'Ovidi o et al.[ 
^004; 



et al. 



Waugh et al. 2006 d'Ovidio et al. 2009 



2008| ). This is mainly due to two reasons: a) tracer ver- 
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tical displacement is usually very small with respect to the hor- 
izontal one; and b) satellite data of any quantity (temperature, 
chlorophyll, altimetry for velocity, etc..) are only available from 
the observation of the ocean surface. 
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Oceanic flows can be considered mainly two-dimensional, 
because there is a great disparity between the horizontal and 
vertical length scales, and they are strongly stratified due to 
the Earth's rotation. There are, however, areas in the ocean 
where vertical motions are fundamental. Firstly there are the 
so-called up welling regions, which are the most biologically 
jpctive marin e zones in the world (Rossi et al. 2008[|Pauly and| 
[Christensen 1995 ). The reason is that due to an Ekmann pump- 
ing mechanism close to the coast, there is a surface uprising of 
deep cold waters rich in nutrients, inducing a high proliferation 
of plankton concentration. Typically, vertical velocities in up- 
welling regions are much larger than in open ocean, but still 
one order of magnitude smaller than horizontal velocities. An- 
other example where there are significant vertical processes are 
mesoscale eddies producing submesoscale structures (fronto- 
genesis), which are responsible for strong ageostrophic vertical 
process, in addition to the vertical exchange thought to occur 
at the eddy interior (Klein and Lapeyre 2009). Thus, the iden- 
tification of the three-dimensional (3d) LCSs in these areas is 
crucial, as well as understanding their correlations with biolog- 
ical activity. Another reason to include the third dimension in 
LCS studies is to investigate the vertical variation in their prop- 
erties. 

The main objective of this paper is the characterization of 
3d LCSs, extracted in an upwelling region, the Benguela area 
in the Southern Atlantic Ocean. For this goal we use Finite- 



83 Size Lyapunov Exponents (FSLEs). FSLEs ( Aurell et al. 1997 



Artale et aL| |1997j measure the separation rate of fluid parti- 
cles between two given distance thresholds. LCSs are com- 
puted as the ridges of the FSLE field (d'Ovidio et al., 2004', 
Molcard et aL| |2006l [Haza et al.||2008, d'Ovidio et al , 2009; 



Poje et al.|[20TQl|Haza et al.|[20TQl ). The rigorous definition of 
LCS as ridges of a Lagrangian stretching measure was given for 
the Finite-Time Lyapunov Exponents (FTLE) in Shadden et al. 



(2005) and Lekien et al. (2007), which are closely related to 
FSLEs. More recently, hyperbolic LCS have been defined in- 
dependently of such stretching measures by Haller p011|). Fol- 
owing many previous studies ( d'Ovidio et aH|2004{" Molcard 



[etaL] [20061 [d'Ovidio al.||2009HBranicki and Wiggins ||2009| ) 
we adopt the mathematical results for Finite-Time Lyapunov 
Exponents (FTLE) to FSLE, assuming them to be valid. In par- 



ticular, we assume that LCS are identified with ridges ( Haller 



2001), i.e., the local extrema of the FTLE field, and also we ex- 



pect, in accordance to the results in [Shadden et al.[ p005[ ) and 
Lekien et"aL] ^2007*) for FTLEs, that the material flux through 



these LCS is small and that they are transported by the flow as 
quasi-material surfaces. 

To confirm that our identification of LCSs with ridges of the 
FSLE field, we perform (in Sect. Ill) direct particle trajectory 
integrations that show that the computed LCS really organize 
the tracer flow. In our work, we will emphasize the numerical 
methodology since up to now FSLEs have only been compute d 
for the marine surface (an exception is Qzgokmen et al. ( 2011 )). 
We then focus on a particular eddy very prominent in the area 
at the chosen temporal window and study the stirring and mix- 
ing on it's vicinity. Some previous results for Lagrangian ed- 
[jiies were obtained by [Branicki and Kirwan[ ( [2010| ) and [BranickT 



et al.[ p011[ ), applying the methodology of lobe dynamics and 
the turnstile mechanism to eddies pinched off from the Loop 
Current. In this paper we focus on FSLE fields and the asso- 
ciated particle trajectories to study transport in and out of the 
chosen mesoscale eddy. Since this is a first attempt to study 
3d oceanic LCS, more general results (on Benguela and other 
upwelling regions) are left for future work. 

To circumvent the lack of appropriate observational data in 
the vertical direction, we use velocity fields from a numeri- 
cal simulation. They are high resolution simulations from the 
ROMS model (see section |2] below) thus appropriate to study 
regional-medium scale basins. 

The paper is organized as follows: In section II we describe 
the data and methods. In section III we present our results. 
Section IV contains a discussion of the results and Section V 
summarizes our conclusions. 

2. Data and Methods. 

2.7. Velocity data set. 

The Benguela ocean region is situated off' the west coast 
of southern Africa. It is characterized by a vigorous coastal 
upwelling regime forced by equatorward winds, a substantial 
mesoscale activity of the upwelling front in the form of eddies 
and filaments, and also by the northward drift of Agulhas ed- 
dies. 

The velocity data set comes from a regional ocean model 
simulation of the Benguela Region ( Le Vu et aTl[2011[ ). ROMS 
( [Shchepetkin and McWilhamsl [2003, 2005] ) is a split-explicit 
free-surface, topography following model. It solves the incom- 
pressible primitive equations using the Boussinesq and hydro- 
static approximations. Potential temperature and salinity trans- 
port are included by coupling advection/diff'usion schemes for 
these variables. The model was forced with climatological data. 
The data set area extends from 12°S to 35°S and from 4°E to 
19°E (see Fig. 1). The velocity field u = (w, v, w) consists of 
two years of daily averaged zonal (u), meridional (v), and verti- 
cal velocity (w) components, stored in a three-dimensional grid 
with an horizontal resolution of 1/12 degrees ~ 8 km, and 32 
vertical terrain-following levels using a stretched vertical co- 
ordinate where the layer thickness varies, increasing from the 
surface to the ocean interior. Since the ROMS model considers 



ij:he hydro static approximation it is important to note that [Ma- 
piadevan (2006), when comparing results from non-hydrostatic 
and hydrostatic versions of the same model of vertical motions 
at submesoscale fronts, found that while instantaneous vertical 
velocities structures diff'er, the averaged vertical ffux is similar 
in both hydrostatic and non-hydrostatic simulations. 

2.2. Finite-Size Lyapunov Exponents. 

In order to study non-asymptotic dispersion processes such 
as stretching at finite scales and time intervals, the Finite Size 
Lyapunov Exponent (Aurell et al. 1997[ Artale et al. 1997[ ) is 
particularly well suited. It is defined as: 



1 Sf 
^= -log/, 

T ^0 



(1) 



166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 



190 
191 
192 
193 
194 
195 
196 




35°S-^ 



Figure 1: Benguela ocean region. The velocity field domain is limited by the 
continuous black line. The FSLE calculation area is limited by the dash-dot 
black line. Bathymetric contour lines are from ETOPOl global relief model 21 6 
fAmante and EakinsH2009jf starting a m depth up to 4000 m at 500 m interval. 217 



where r is the time it takes for the separation between two par- 
ticles, initially ^0, to reach df. In addition to the dependence 
on the values of ^0 and the FSLE depends also on the initial 
position of the particles and on the time of deployment. Lo- 
cations (i.e. initial positions) leading to high values of this 
Lyapunov field identify regions of strong separation between 
particles, i.e., regions that will exhibit strong stretching during 
evolution, that can be identified with the LCS dBofFetta et al.l 



2Q0T]|d'Qvidio et al.||2QQ4| [Joseph and Legras|[2002l ). 



In principle, for computing FSLEs in three dimensions one 



just needs to extend the method of |d'Qvidio et al.|(2004 ) , that is, 
one needs to compute the time that fluid particles initially sep- 
arated by ^0 = [(^-^0)^ + (^Jo)^ + (feo)^]^^^ need to reach a final 
distance of ^/ = [{dXff + {6yff + {6zffY'^. The main difficulty 
in doing this is that in the ocean vertical displacements (even in 
up welling regions) are much smaller than the horizontal ones, 
and so do not contribute significantly to total particle dispersion 
( Ozgokmen et al. 201 1 ). By the time the horizontal particle dis- 



persion has scales of tenths or hundreds of kilometers (typical 
mesoscale structures are studied using df ^ lOO^m (d'Ovidio 



et al.[|2004] )), particle dispersion in the vertical can have at most 



scales of hundreds of meters and usually less. This means that 
the vertical separation will not contribute significantly to the ac- 
cumulated distance between particles. In addition, since length 
scales in the horizontal and vertical diff'er by several orders of 
magnitude, one faces the impossibility of assigning equal ^0 
to the horizontal and vertical particle pairs. It should be noted 
however that these shortcomings arise from the diff'erent scales 
of length and time that characterize horizontal and vertical dis- 
persion processes in the ocean, and so should not be seem as 
intrinsic limitations of the method. For non-oceanic flows a di- 



rect generalization of FSLEs is straightforward. 

Thus, in this paper we implemented a quasi three- 
dimensional computation of FSLEs. That is, we make the com- 
putation for every (2d) ocean layer, but where the particle tra- 
jectories calculation use the full 3d velocity field. I.e., at each 
level (depth) we set dzo = 0, and the final distance is computed 
without taking the vertical distance between particles. It is im- 
portant to note that, since we allow the particles to evolve in the 
full 3d velocity field, we take into account vertical quantities 
such as vertical velocity shear that may influence the horizontal 
separation between particle pairs. 

There are other possible approaches to the issue of difl'er- 
ent scales in the vertical and horizontal. One way is to assign 
anisotropic initial and final displacements in the FSLE calcula- 
tion (i. e., including a Szo and Szf much smaller than the hor- 
izontal initial and final separations). A second approach is to 
use diff'erent weights for the horizontal and vertical separations 
in the calculations of the distance, perhaps in combination with 
the first. We have cheked both alternatives and found that, with 
reasonable choices of initial and final distances and distance 
metrics, the results were equivalent to the quasi-3d computa- 
tion. The reason is that actual dispersion is primarily horizontal 
as commented above. 



More in detail, a grid of initial locations xq in the longi- 
tude/latitude/depth geographical space (0, 6, z), fixing the spa- 
tial resolution of the FSLE field, is set up at time t. The horizon- 
tal distance among the grid points, ^0, was set to 1/36 degrees 
3 km), i.e. three times finer resolution than the velocity field 
( [Hernandez-Carrasco et al.[ |2Q11| ), and the vertical resolution 
(distance between layers) was set to 20 m in order to have a 
good representation of the vertical variations in the FSLE field. 
Particles are released from each grid point and their three di- 
mensional trajectories calculated. The distances of each par- 
ticle with respect to the ones that were initially neighbors at 
an horizontal distance ^0 are monitored until one of the hori- 
zontal separations reaches a value By integrating the three 
dimensional particle trajectories backward and forward in time, 
we obtain the two diff'erent types of FSLE maps: the attract- 
ing LCS (for the backward), and the repelling LCS (forward) 



( d'Ovidio et al.| [2004, Joseph and Legras , 2002) . We obtain in 
this way FSLE fields with a horizontal spatial resolution given 
by ^0- The final distance Sf was set to 100 km, which is, as 
already mentioned, a typical length scale for mesoscale studies. 
The trajectories were integrated for a maximum of T = 178 
days (approximately six months) using an integration time step 
of 6 hours. When a particle reached the coast or left the veloc- 
ity field domain, the FSLE value at its initial position and initial 
time was set to zero. If the interparticle horizontal separation 
remains smaller than 5f during all the integration time, then the 
FSLE for that location is also set to zero. 



247 The equations of motion that describe the evolution of parti- 



3 



248 cle trajectories are 
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dt 
do 
dt 
dz 
dt 



1 u((f),0,zj) 
cos(e) ' 

w(0, 0, z, t). 



(2) 
(3) 
(4) 



where is longitude, 6 is latitude and z is the depth. is the 
radial coordinate of the moving particle = R - z, with R = 
6371 km the mean Earth radius. For all practical purposes, R^ ~ 
R. Particle trajectories are integrated using a 4^^ order Runge- 
Kutta method. For the calculations, one needs the (3d) velocity 
values at the current location of the particle. Since the six grid 
nodes surrounding the particle do not form a regular cube, direct 
trilinear interpolation can not be used. Thus, an isoparametric 
element formulation is used to map the nodes of the velocity 
grid surrounding the particles position to a regular cube, and 



This is done in the next section by comparing actual particle tra- 
jectories with the computed LCSs. 

Differently than 2d, where LCS can be visually identified as 
the maxima of the FSLE field, in 3d the ridges are hidden within 
the volume data. Thus, one needs to explicitly compute and 
extract them, using the definition of LCSs as the ridges of the 
FSLEs. A ridge L is a co-dimension 1 orientable, diff'erentiable 
manifold (which means that for a three-dimensional domain D, 



262 an inverse isoparametric mapping scheme ( Yuan et al. 1994) 



263 is used to find the coordinates of the interpolation point in the 

264 regular cube coordinate system. 

265 2.3. Lagrangian Coherent Structures. 



In 2d, LCS practically coincide with (finite-time) stable and 321 
unstable manifolds of relevant hyperbolic structures in the flow 322 



(Haller| [2Q00l [Hafler and Yu^ [2QQQ1 [Joseph and Legras 



2002]k The structure of these last objects in 3d is generally 



jidge s are surfaces) satisfying the following conditions ( [Lekien 
|et^|2007j : 

1 . The field A attains a local extremum at L. 

2. The direction perpendicular to the ridge is the direction of 
fastest descent of A at L. 

Mathematically, the two previous requirements can be ex- 
pressed as 

n'^V^ = 0, (5) 

n^Hn = min u^Hu < 0, (6) 

llulNl 

where VA is the gradient of the FSLE field A, n is the unit nor- 
mal vector to L and H is the Hessian matrix of A. 

The method used to extract the ridges from the scalar field 
^(xq, t) is from Schultz et al. (2010). It uses an earlier ( |Eberly| 
[et al. 1994) definition of ridge in the context of image analy- 
sis, as a generalized local maxima of scalar fields. For a scalar 
field / : ^ R with gradient g = V/ and Hessian H, a 
J-dimensional height ridge is given by the conditions 



much more complex than in 2d ( |Haller|[2001[ Pouransari et al. 
'2010), and they can be locally either lines or surfaces. As com- 325 
mented before, however, vertical motions in the ocean are slow. 326 
Thus, at each fluid parcel the strongest attracting and repelling 327 
directions should be nearly horizontal. This, combined with 
the incompressibility property, implies that the most attracting 
and repelling regions (i.e. the LCSs) should appear as almost 329 
vertical surfaces, since the attraction or repulsion should occur 330 
normally to the LCS. As a consequence, the LCSs will have a 331 
"curtain-like" geometry, with deviations from the vertical due 332 
to either the orientation of the most attracting or repelling direc- 333 
tion deviating from the horizontal, or when strong vertical shear 334 
produces variations along the vertical in the most repelling or 335 
attracting regions in the flow. We expect the LCS sheet-like ob- 336 
jects to coincide with the strongest hyperbolic manifolds when 337 
these are two dimensional, and to contain the strongest hyper- 333 
bolic lines. 339 
The curtain-like geometry of the LCS was already com- 340 
jnented in [Branicki and Malek-Madani| ( [2010| ), [Branicki and] 341 
|Kirwan| ([2010), or 'B ranicki et al. ([2011). In the latter paper 342 
it was shown that, in a 3d flow, these structures would appear 343 
mostly vertical when the ratio of vertical shear of the horizon- 344 
tal velocity components to the average horizontal velocities is 345 
small. This ratio also determines the vertical extension of the 346 
structures. In Branicki and Kirwan ( 2010| ), the argument was 347 
used to construct a 3d picture of hyperbolic structures from the 348 
computation in a 2d slice. In the present paper we confirm the 349 
curtain-like geometry of the LCSs, and show that they are rele- 350 
vant to organize the fluid flow in this realistic 3d oceanic setting. 351 



T 



and ai < 0, 



(7) 



where aiJ e {1,2,..., n}, are the eigenvalues of H, ordered 
such that ai > . . . > an, and C/ is the eigenvector of H associ- 
ated with ai. For n = becomes 

g'^es = and 0^3 < 0. (8) 

This ridge definition is equivalent to the one given by ([5]) since 
the unit normal n is the eigenvector (when normalized) associ- 
ated with the minimum eigenvalue of H. In other words, in R^ 
the ei , 62 eigenvectors point locally along the ridge and the 63 
eigenvector is orthogonal to it. 

The ridges extracted from the backward FSLE map approxi- 
mate the attracting LCS, and the ridges extracted from the for- 
ward FSLE map approximate the repelling LCS. The attract- 
ing ones are the more interesting from a physical point of view 
( d'Qvidio et al.[ |2004[ |2009| ), since particles (or any passive 
scalar driven by the flow) typically approach them and spread 
along them, giving rise to filament formation. In the extrac- 
tion process it is necessary to specify a threshold s for the ridge 
strength {asl so that ridge points whose value of as is lower 
(in absolute value) than s are discarded from the extraction pro- 
cess. Since the ridges are constructed by triangulations of the 
set of extracted ridge points, the s threshold greatly determines 
the size and shape of the extracted ridge, by filtering out re- 
gions of the ridge that have low strength. The reader is referred 
to Schultz et al. ( 2010| for details about the ridge extraction 



method. The height ridge definition has been used to extract 
LCS from FTLE fields in several works (see, among others, 
ISadlo and Pei"totl ( |2007| ). 
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Figure 2: Vertical profile of 30 day average backward and forward FSLE. 
The 30 day average field was spatially averaged at each layer over the FSLE 
calculation area to produce the vertical profiles. The backward FSLE average 
is shown in continuous and the forward FSLE is shown in dashed. 

352 3. Results 

382 

353 3.1. Three dimensional FSLE field 

354 The three dimensional FSLE field was calculated for a 30 384 

355 day period starting September 17, with snapshots taken every 385 

356 2 days. The fields were calculated for an area of the Benguela 386 

357 ocean region between latitudes 20°S and 30°S and longitudes 387 

358 8°E to 16°E (see figure [T]). The area is bounded at NW by the 388 

359 Walvis Ridge and the continental slope approximately bisects 389 

360 the region from NW to SE. The western half of the domain has 390 

361 abyssal depths of about 4000 m. The calculation domain ex- 391 

362 tended vertically from 20 up to 580 m of depth. Both backward 392 

363 and forward calculations were made in order to extract the at- 393 

364 tracting and repelling LCS. 394 

365 Figure |2] displays the vertical profile of the average FSLE 395 

366 for the 30 day period. There are small diff'erences between the 396 

367 backward and the forward values due to the diff'erent intervals 397 

368 of time involved in their calculation. But both profiles have a 393 

369 similar shape and show a general decrease with depth. There 399 

370 is a notable peak in the profiles at about 100 m depth that indi- 400 

371 cates increased mesoscale variability (and transport, as shown 

372 in Sect. [3^at that depth). 

373 A snapshot of the attracting LCSs for day 1 of the calculation 

374 period is shown in figure |3] As expected, the structures appear 402 

375 as thin vertical curtains, most of them extending throughout the 403 

376 depth of the calculation domain. The area is populated with 404 

377 LCS, denoting the intense mesoscale activity in the Benguela 405 

378 region. As already mentioned, in three dimensions the ridges 406 

379 are not easily seen, since they are hidden in the volume data. 407 

380 However the horizontal slices of the field in figure [3] show that 408 

381 the attracting LCS fall on the maximum backward FSLE field 409 
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Figure 3: Attracting LCS (blue) for day 1 of the calculation period, together 
with horizontal slices of the backward FSLE field at 120 m and 300 m depth. 
Colorbar refers to colormap of horizontal slices. The units of the colorbar are 

lines of the 2d slices. The repelling LCS (not shown) also fall 
on the maximum forward FSLE field lines of the 2d slices. 

Since the A value of a point on the ridge and the ridges 
strength 0^3 are only related through the expressions ^ and ([8]), 
the relationship between the two quantities is not direct. This 
creates a difficulty in choosing the appropriate strength thresh- 
old for the extraction process. A too small value of s will result 
in very small LCS that appear to have little influence on the 
dynamics, while a greater value will result in only a partial ren- 
dering of the LCS, limiting the possibility of observing their 
real impact on the flow. Computations with several values of s 
lead us to the optimum choice s = 2Q day~^m~^, meaning that 
grid nodes with (23 < -20 day~^m~^ were filtered out from the 
LCS triangulation. 

We have seen in this section an example of how the ridges 
of the 3d FSLE field, the LCS, distribute in the Benguela ocean 
region. Their ubiquity shows their impact on the transport and 
mixing properties. In the next section we concentrate on the 
properties of a single 3d mesoscale eddy. 

3.2. Study of the dynamics of a relevant mesoscale eddy 

Let us study a prominent cyclonic eddy observed in the data 
set. The trajectory of the center of the eddy was tracked and it 
is shown in figure |4] The eddy was apparently pinched off at 
the upwelling front. At day 1 of the FSLE calculation period 
its center was located at latitude 24.8°S and longitude 10.6°E, 
leaving the continental slope, and having a diameter of approx- 
imately 100 km. One may ask: what is its vertical size? is it 
really a barrier, at any depth, for particle transport? 
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particle transport: eddies can be defined as the region bounded 
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Figure 4: Trajectory (advancing from NE to SW) of the eddy center inside the ^eo 
calculation domain. Circles indicate the center location during the 30 day FSLE 46i 
calculation period, and squares previous and posterior positions. Bathymetric 
lines same as in figurefTl 
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To properly answer these questions the eddy, in particular its 
frontiers, should be located. From the Eulerian point of view 
it is commonly accepted that eddies are delimited by closed 
contours of vorticity and that the existence of strong vorticity 
gradients prevent the transport in and out of the eddy. Such . 
transport may occur when the eddy is destroyed or undergoes 
strong interactions with other eddies ( Provenzale 1999| ). In a 
Lagrangian view point, however, an eddy can be defined as a re- 
gion delimited by intersections and tangencies of LCS, whether 
jn 2d or 3d space. The eddy itself is an elliptic structure ( Haller| 
[and YuanI |25Q01 [Branicki and Kirwanl [MTOl [Branicki et al 



2011 ). In this Lagrangian view of an eddy, the transport inhi- 



bition to and from the eddy is now related to the existence of 
these transport barriers delimiting the eddy region, which are 
known to be quasi impermeable. 

Using the first approach, i.e., the Eulerian view, the vertical 
distribution of the Q-criteria (Hunt e t aL[|1988[|Jeo ng and Hus-| 
sain , 1995 ) was used to determine the vertical extension of the 
mesoscale eddy. The Q criterium is a 3d version of the Okubo- 
Weiss criterium (Qkubo| [T970l |Weiss| [T99T] ) and measures the 
relative strength of vorticity and straining. In this context, ed- 
dies are defined as regions with positive Q, with Q the second 
invariant of the velocity gradient tensor 



e = i(lin|p-lis|p). 



intersecting or tangent repelling and attracting LCS ( |Bran- 
icki and Kirwan[ |2010[ Branicki et al. 201 1| ). Using this cri- 



terion, and first looking at the surface located at 200 m depth, 
we see in Fig. |6]that certainly the Eulerian eddy seems to be 
located inside the area defined by several intersections and tan- 
gencies of the LCS. This eddy has an approximate diameter of 
100 km. In the south-north direction there are two intersections 
that appear to be hyperbolic points (HI and H2 in figure [6]). 
In the West-East direction, the eddy is closed by a tangency at 
the western boundary, and a intersection of lines at the eastern 
boundary. The eddy core is devoid of high FSLE lines, indicat- 
ing that weak stirring occurs inside (d'Qvidio et al. 2004| ). As 
additional Eulerian properties, we note that near or at the inter- 
sections HI and H2 the Q-criterium indicates straining motions. 
In the case of H2, figure [5] (right panel) indicates high shear up 
to 200 m depth. The fact that the hyperbolic regions HI and 
H2 lie in strain dominated regions of the flow (2 < 0) high- 
lights the connection between hyperbolic particle behavior and 
instantaneous hyperbolic regions of the flow. The ridges of the 
FSLE field, however, do not remain in the negative Q regions 
but cross into rotation dominated regions with 2 > 0. This indi- 
cates that there are some diff'erences between the Eulerian view 
(Q) and the Lagrangian view (FSLE). It is the latter that can 
be understood in terms of particle behaviour as limiting regions 
of initial conditions (particles) that stay away from hyperbolic 



regions for long enough time (Halle r and Yuan|[200 ). 

In 3d, the eddy is also surrounded by a set of attracting and 
repelling LCS (figure [7]), calculated as explained in Subsection 
|2.3| The lines identified in figure [6] are now seen to belong to 
the vertical of these surfaces. 

Note that the vertical extent of these surfaces is in part de- 
termined by the strength parameter used in the LCS extraction 
process, so their true vertical extension is not clear from the 
results presented here. On the south, the closure of the La- 
grangian eddy boundary extends down to the maximum depth 
of the calculation domain, but moving northward it is seen that 
the LCS shorten their depth. Probably this does not mean that 
the eddy is shallower in the North, but rather that the LCS are 
losing strength (lower \a2,\) and portions of it are filtered out 
by the extraction process. In any case, it is seen that as in 
two-dimensional calculations, the LCS delimiting the eddy do 



jiot per fectly coincide with its Eulerian boundary ( [Joseph and 
[Legrasl |2002| ), and we expect the Lagrangian view to be more 
relevant to address transport questions. 

In the next paragraphs we analyze the fluid transport across 
the eddy boundary. Some previous results for Lagrangian ed- 
ies were obtained by Branicki and Kirwan (2010') and Branicki 



(9) 489 

490 

where = tri^l^^), ||S|p = triSS^) and O, S are the anti- 491 
symmetric and symmetric components of Vu. Using 2 = as 492 
the Eulerian eddy boundary, it can be seen from Fig. [5] that the 493 
eddy extends vertically down to, at least, 600 m. 494 
Let us move to the Lagrangian description of eddies, which 495 
is much in the spirit of our study, and will allow us to study 496 



l et al.|(2011) . Applying the methodology of lobe dynamics and 
the turnstile mechanism to eddies pinched ofl' from the Loop 
Current, 'Bran icki and Kirwan] ( [201 0| observed a net fluid en- 
trainment near the base of the eddy, and net detrainment near 
the surface, being fluid transport in and out of the eddy essen- 
tially confined to the boundary region. Let us see what happens 
in our setting. 

We consider six sets of 1000 particles each, that were re- 
leased at day 1 of the FSLE calculation period, and their trajec- 
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Figure 5: Colormap of g-criterium. White contours have Q = 0. Day 1 of the 30 day FSLE calculation period. Left panel: Latitude 24.5° S ; Rigth panel: Longitude 
10.5° E. Colorbar values are Q x 10^^ s'^. 




Figure 6: Q-criterium map at 200 m depth together with patches of backward 
(blue) and forward (green) FSLE values. Black dashed lines have Q = 0. FSLE 
patches contain the highest 60% of FSLE values. Colorbar values are Q x 
]^q10 ^-2 jjig Q^^y study is the clear region in between points HI and H2. 



Figure 7: 3d LCSs around the mesoscale eddy at day 1 of the 30 day FSLE 
calculation period. Green: repelling LCS; Blue: attracting LCS. 
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497 tories integrated by a fourth-order Runge-Kutta method with a 521 

498 integration time step of 6 hours. The sets of particles were re- 522 

499 leased at depths of 50, 100, 200, 300, 400 and 500 m. In figure 523 

500 [8] we plot the particle sets together with the Lagrangian bound- 524 

501 aries of the mesoscale eddy viewed in 3d. A top view is shown 525 

502 in figure [9] As expected, vertical displacements are small. 526 

503 At day 3 (top left panel of figures [8] and [9]) it can be seen that 527 

504 there is a diff'erential rotation (generally cyclonic, i.e. clock- 528 

505 wise) between the sets of particles at diff'erent depths. The shal- 529 

506 lower sets rotate faster than the deeper ones. This diff'erential 530 

507 rotation of the fluid particles could be viewed, in a Lagrangian 531 

508 perspective, as the fact that the attracting and repelling strength 532 

509 of the LCS that limit the eddy varies with depth. Note that the 533 

510 six sets of particles are released at the same time and at the same 534 

511 horizontal positions, and thereby their difl'erent behavior is due 535 

512 to the variations of the LCS properties along depth. sse 

513 At day 13 the vortex starts to expel material trough filamen- 537 

514 tation (Figs [8] and |9] top right panels). A fraction of the par- 533 

515 tides approach the southern boundaries of the eddy from the 539 

516 northeast. Those to the west of the repelling LCS (green) turn 540 

517 west and recirculate inside the eddy along the southern attract- 541 

518 ing LCS (blue). Particles to the east of the repelling LCS turn 542 

519 east and leave the eddy forming a filament aligned with an at- 543 

520 trading (blue) LCS. At longer times trajectories in the south of 544 



the eddy are influenced by additional structures associated to 
a difl'erent southern eddy. At day 29 (bottom right panels) the 
same process is seen to have occurred in the northern boundary, 
with a filament of particles leaving the eddy along the northern 
attracting (blue) LCS. The filamentation seems to begin earlier 
at shallower waters than at deeper ones since the length of the 
expelled filament diminishes with depth. However all of the 
expelled filaments follow the same attracting LCS. Figure 10 
shows the stages previous to filamentation in which the LCS 
structure, their tangencies and crossings, and the paths of the 
particle patches are more clearly seen. Note that the LCS do not 
form fully closed structures and the particles escape the eddy 
through their openings. The images suggest lobe-dynamics pro- 
cesses, but much higher precision in the LCS extraction would 
be needed to really see such details. 

This filamentation event seems to be the only responsible for 
transport of material outside of the eddy, since the rest of the 
particles remained inside the eddy boundaries. To get a rough 
estimate of the amount of matter expelled in the filamentation 
process we tracked the percentage of particles leaving a circle 
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of diameter 200 km centered on the eddy center. In Fig. 
the time evolution of this percentage is shown for the particle 
sets released at diff"erent depths. The onset of filamentation is 
clearly visible around days 9-12 as a sudden increase in the per- 
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Figure 9: Top view of the evolution of particle patches and LCSs shown in Fig. 8. Top left: day 3; Top right: day 13; Bottom left: day 19: Bottom right: day 29. 
Colors as in figure 8. 



545 centage of particles leaving the eddy. The percentage is maxi- 566 

546 mum for the particles located at 100 m depth and decreases as se? 

547 the depth increases. At 400 and 500 m depth there are no parti- 568 

548 cles leaving the circle. There is a clear lag between the onset of 569 

549 filamentation between the different depths: the onset is simulta- 570 

550 neous for the 40 m and 100 m depths but occurs later for larger 571 

551 depths. 572 

573 

552 4. Discussion. 

575 

553 The spatial average of FSLEs defines a measure of stirring 576 

554 and thus of horizontal mixing between the scales used for its 577 

555 computation. The larger the average, the larger the mixing ac- 573 

2004| ). The general trend in the verti- 



556 tivity (d'Ovidio et al. 



cal profiles of the average FSLE (Fig. [3]) shows a reduction of 530 
mesoscale mixing with depth. There is however a rather inter- 531 
esting peak in this average profile occurring at 100 m, i.e. close 532 
to the thermocline. It could be related to submesoscale pro- 533 
cesses that occur alongside the mesoscale ones. Submesoscale 534 
is associated to filamentation (the thickness of filaments is of 535 
the order of 10 km or less), and we have seen that the filamen- 536 



tation and the associated transport intensity (Fig. 1 1 ) is higher 537 
at 100 m depth. It is not clear at the moment what is the precise 533 



mechanism responsible for this increased activity at around 100 
m depth (perhaps associated to instabilities in the mixed layer), 
but we note that the intensity of shearing motions (see the Q 
plots in|5]) is higher in the top 200 meters. Less intense filamen- 
tation could be caused by reduction of shear in depths larger 
than these values. 

From an Eulerian perspective, it is thought that vortex fil- 
amentation occurs when the potential vorticity (PV) gradient 
aligns itself with the compressional axis of the velocity field, 
in strain coordinates ( Louazel and Hua ( 2004| ): Lapeyre et al. 
( [1999] )). This alignment is accompanied by exponential growth 
of the PV gradient magnitude. The fact that the filamentation 
occurs along the attracting LCS seems to indicate that this ex- 
ponential growth of the PV gradient magnitude occurs across 
the attracting LCS. 

In the specific spatiotemporal area we have studied, and in 
particular, for the eddy on which we focus sed our analysis, we 
have confirmed that the structure of the LCSs is "curtain-like", 
so that the strongest attracting and repelling structures are qua- 
sivertical surfaces. Their vertical extension would depend of the 
physical transport properties, but it is also altered by the partic- 
ular threshold parameter selected to extract the LCSs. These 
observations imply that transport and stirring occurs mainly on 
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Figure 10: Top view of the initial stages of evolution of the particle patches and LCSs of Figs. 8 and 9. Top left: day 7; Top right: day 9; Bottom left: day 11: 
Bottom right: day 13. Colors as in figure 8. 
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the horizontal, which is a reasonable result considering the dis- 612 
parity between horizontal and vertical velocities in the ocean, 613 
and its stratification. However, we should mention that our re- 6i4 
suits are not fully generalizable to all ocean situations, and that 615 
any ocean area or oceanic event should be studied in particular 6I6 
to reveal the shape of the associated 3d LCS. 617 
Some comments follow about the nature of vertical trans- 6I8 
port structures. FSLEs are suited to the identification of hyper- 619 
bolic structures (structures that exhibit high rates of transversal 620 
stretching or compression in their vicinity). The question is 621 
if one can expect that structures responsible for vertical trans- 622 
port will also exhibit substantial (vertical) stretching. This is 623 
not so clear in the ocean for the reasons already indicated. If 624 
one considers the case (relevant to our work) of purely isopy- 625 
cnal flow, then strong vertical stretching would be associated 
with a rapid divergence of isopy cnic surfaces. In the case of 
coastal upwelling, for instance, the lifted isopycnic surfaces 
move vertically in a coherent fashion, so one should not expect 627 
strong vertical divergence of particles flowing along neighbour- 628 
ing isopycnic surfaces. This is just an example of the fact that it 629 
is possible that coherent vertical motions do not imply the pres- 63o 
ence of hyperbolic coherent structures such as those the FSLE 63i 
may indicate. 632 



Another possible limitation worth mentioning is the velocity 
field resolution and its relation to the intensity of the vertical 
velocity. It is accepted that in fronts or in the eddy periphery, 
vertical velocities are significantly greater than, for instance, in 
the eddy interior. These zones of enhanced vertical transport 
correspond to submesoscale features that were not adequately 
captured in the velocity field used in this work due to its coarse 
resolution, since submesoscale studies usually have resolutions 
< 10 km (the literature on this subject is quite large, so we refer 
the reader to Klein and Lapeyre| ( |2009| ) and Levy ( 2008| ) ). 

In any case, a most important point for the LCS we have 
computed is that in 3d, as in 2d, they act as pathways and bar- 
riers to transport, so that they provide a skeleton organizing the 
transport processes. 

5. Conclusions 

Three dimensional Lagrangian Coherent Structures were 
used to study stirring processes leading to dispersion and mix- 
ing at the mesoscale in the Benguela ocean region. We have 
computed 3d Finite Size Lyapunov Exponent fields, and LCSs 
were identified with the ridges these fields. LCSs appear as 
quasivertical surfaces, so that horizontal cuts of the FSLE fields 
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665 mesoscale transport and mixing processes. 
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Figure 1 1 : Percentage of particles outside a 200 km diameter circle centered at 
the eddy center, as a function of time. 



gives already a quite accurate vision of the 3d FSLE distribu- 
tion. These quasivertical surfaces appear to be coincident with 
the maximal lines of the FSLE field (see fig. [3]) so that sur- 
face FSLE maps could be indicative of the position of 3d LCS, 
as long as the vertical shear of the velocity does not result in a 
significant deviation of the LCS with respect to the vertical. Av- 
erage FSLE values generally decrease with depth, but we find 
a local maximum, and thus enhanced stretching and dispersion, 
at about 100 m depth. 

We have also analyzed a prominent cyclonic eddy, pinched 
off the upwelling front and study the filamentation dynamics 
in 3d. Lagrangian boundaries of the eddy were made of in- 
tersections and tangencies of attracting and repelling LCS that 
apparently emanating from two hyperbolic locations North and 
South of the eddy. The LCS are seen to provide pathways and 
barriers organizing the transport processes and geometry. This 
pattern extends down up to the maximum depth were we cal- 
culated the FSLE fields (~ 600 m), but the exact shape of the 
boundary is difficult to determine due to the decrease in ridge 
strength with depth. This caused some parts of the LCS not to 
be extracted. The inclusion of a variable strength parameter in 
the extraction process is an important step to be included in the 
future. 

The filamentation dynamics, and thus the transport out of the 
eddy, showed time lags with increasing depth. This arises from 
the vertical variation of the flow field. However the filamenta- 
tion occurred along all depths, indicating that in reality vertical 
sheets of material are expelled from these eddies. 

Many more additional studies are needed to further clarify 
the details of the geometry of the LCSs, their relationships with 
finite-time hyperbolic manifolds and three dimensional lobe dy- 
namics, and specially their interplay with mesoscale and sub- 
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